function y = funlab5ej4()
	
	N = 100;	
    a = 0;
    b =  pi/2;
    
	y = simpson_compuesta('f',a, b, N)
	
endfunction


function res = simpson_compuesta(fun,a, b, N)

	h = (b-a)/N;
    sum1 = 0;
    sum2 = 0;

    for j = 1:(N/2)-1
		i = 2 * j;
        x = a + (i * h);
        sum1 += feval(fun,x);
    endfor

    for j = 1 : N/2
		i = (2 * j) - 1;
        x = a + (i * h);
        sum2 += feval(fun,x);
    endfor

    res = (h/3) * (feval(fun,a) + 2*sum1 + 4*sum2 + feval(fun,b));
    
endfunction



